1/7 



IR 



110 



CHANNEL 



IR * 




120 



TIME 



VARIABLE 
SIRF 



TIME 



FIG. 1 



200 

S 

RECEIVER MODEM 
230 



110 



1 



COMM. 
PORT(S) 



120 



VARIABLE 
SIRF 




240 



DECODER 



DATA STORAGE DEVICE 



FILTER COEFFICIENT 
DETERMINATION PROCESS 



FIG. 2 



2/7 



300 



FILTER 

COEFFICIENT DETERMINATION 
PROCESS 




310 



INITIALIZE THE SIRF FILTER 
TO AN ARBITRARY VALUE 



INPUT: CHANNEL 
IMPULSE RESPONSE 



320 



DEFINE SETS C 2 , C 4 , C 5 (NON-LINEAR 

FILTER) OR C p C 3 , C 5 (LINEAR FILTER) 

TO SPECIFY DESIRED CHARACTERISTICS 
IN TIME AND FREQUENCY DOMAIN 



350 



DEFINE PROJECTION OPERATORS 
P 2 , P 4 , P 5 OR P 1t P 3 , P 5 



USE PROJECTION OPERATORS P 2 , P 4 , P 5 
TO PROJECT ONTO SETS C 2 , C 4 , C 5 (NON- 
LINEAR FILTER) OR PROJECTION 
OPERATORS P 1t P 3 , P 5 TO PROJECT 

ONTO SETS C p C 3 , C 5 (LINEAR FILTER) 



370 — 



1 


r 


ITERATIVELY CONTINUE PROJECTIONS 
UNTIL INTERSECTION IS REACHED 




r 


oo n GENERATE N 
oou FILTER COEFFICIENTS 




r 



END 



360 — 



FIG. 3 



3/7 



400 

/ 

/ length = 64; //fast fourier transform length 

delta = .1; %input ('Enter stop band tolerance: f ); 
beta = 0.2; %input ('Enter pass band tolerance: '); 
wp = 8; 
ws = 7; 
I mu = 28; %CP length 

/ load channeljmpulse_resp.dat; %load channel impulse response 
/ N = 19; % sirf filter length 
4 1 0 < ch = channeMmpulse_resp; 
\ mm = max(abs(ch)); 
\ ch = 2047*ch./mm; 

NN = size(channeMmpulse_resp); 
M = NN(1); 
for i = 1:mu 

mv(i) = 1 ; %initialize mv 
end 

i C = conv(mv,abs(ch)); 

V [Y,l] = max(abs(C)); %find max which is beginning of Gl (CP) 



/ % construct the channel impulse response matrix 
for i = 1 :M 
forj = l:N 

if i - j + 1 <= 0 

break; 
else 

H(iJ) = ch(i-j + 1); 
end 
end 
end 
n = 1; 

fori = M + 1:M + N-1 
n = n + 1 ; 
k = 0; 
for j = n:N 

k = k + 1 ; 

H(ij) = ch(M-k+1); 

end 
end 

for i = 1 :mu 

mv(i) = 1; 
end 

C = conv(mv,abs(ch)); 
[Yl] = max(abs(C)); 



FIG. 4A 



4/7 



for i = 1 :mu 
mv(i) = 1; 
440 < end 

C = conv(mv,abs(ch)); 
[Y,l] = max(abs(C)); 



/ % initialize the SIRF filter to an arbitrary value 
' c = [ -0.0593 
0.01047 
-0.062386 
0.02418 
-0.065114 
0.030031 
-0.039083 
0.13789 
-0.10266 
-0.014681 
-0.11224 

00546 FIG. 4B 

-0.12642 
0.18608 
-0.020895 
0.38407 
-0.30117 
i -0.37885 
\ 0.42326]; 



g = fft(clength); 
f = abs(g); 
450 < lamda = 1; 
tol = 50; 

cold = sum(abs(c)); 
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/ %begin iteration 

/ for ii = 1:20 % 20 is the number of iterations 
% projection on the set C2 
for m = wp:length/2 
if ( f(m) > (1 + beta)) 

gg = g(m); 

g(m) = (1 + beta) / f(m)*complex(real(g(m)),imag(g(m))); 
g(m) = gg + lamda*(g(m) - gg); 
end 

if ( f(m) < (1 - beta)) 

gg = g(m); 

g(m) = (1 - beta) / f(m)*complex(real(g(m)),imag(g(m))); 
g(m) = gg+lamda*(g(m) - gg); 
end 
end 

for m = 1 :ws 

if (f(m) > delta) 

gg = g(m); 

g(m) = delta / f(m)*complex(real(g(m)),imag(g(m))); 
g(m) = gg+lamda*(g(m) - gg); 
end 
end 

for m = 2:length/2 

g(m + length/2) = conj(g(length/2 - m + 2)); 
\ end 

\ cr = real(ifft((g),length)); %transform to time domain 
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% time domain projection on C4 set 
for i = 1:N 

c(i) = cr(i); %project into the set C4 

end 

% time domain projection on C5 set 
for n = I - mu - 10:1 - mu 
norm = 0; 
prod = 0; 
%prod1 = 0; 
for i = 1:N 

norm = norm + H(n,i) A 2; 

prod = prod + c(i)*H(n,i); 

%prod1 = prodl + c(i)*ch(n - i + 1); 
end 

if (prod > toi) 
for nn = 1:N 
cc = c(nn); 

c(nn) = c(nn) + ((to! - prod /norm) * H(n,nn); 
c(nn) = cc + lamda*(c(nn) - cc); 
end 
end 

if (prod < -tol) 
for nn = 1:N 
cc = c(nn); 

c(nn) = c(nn) + ((-tol - prod /norm) * H(n,nn); 
c(nn) = cc + lamda*(c(nn) - cc); 
end 
end 



forn = l;l + 100 
norm = 0; 
prod = 0; 
fori = 1:N 

norm = norm + H(n,i) A 2; 

prod = prod + c(i)*H(n,i); 
end 

if (prod > tol) 
for nn = 1:N 
cc = c(nn); 

c(nn) = c(nn) + ((tol - prod) /norm) * H(n,nn); 
c(nn) = cc + lamda*(c(nn) - cc); 

end 
end 

if (prod < -toi) 
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/" fornn = 1:N 

cc = c(nn); 

c(nn) = c(nn) + ((-tol - prod) /norm * H(n,nn); 
c(nn) = cc + Iamda*(c(nn) - cc); 
end 
end 

end 

ss = abs(co!d - sum(abs(c))) 
cold = sum(abs(c)); 

g = fft(cjength); %transform to frequency domain 
f = abs(g); 
y end 

f fori = 1:N 

480 < sirf(i) = c(i); %sirf will hold the SIRF filter coefficients 
end 
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